Subdiffusion-limited reactions 
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We consider the coagulation dynamics A + A — * A and 
A + A ^ A and the annihilation dynamics A + A — ► for par- 
ticles moving subdiffusively in one dimension. This scenario 
combines the "anomalous kinetics" and "anomalous diffusion" 
problems, each of which leads to interesting dynamics sepa- 
rately and to even more interesting dynamics in combination. 
Our analysis is based on the fractional diffusion equation. 

PACS numbers: 02.50.Ey,82.40.-g,05.60.-k,05.70.Ln 



Diffusion-limited reactions in constrained geome- 
tries bave been studied intensely because they exhibit 
"anomalous kinetics," that is, behavior different from 
that predicted by the laws of mass action in well-stirred 
systems [jj]. Among the simplest and most extensively 
studied are single species diffusion-limited coagulation 
(A + A ^ A or A + A ^ A) || and annihilation 
(A + A — ► 0) (|||. These reactions, which show anoma- 
lous behavior in one dimension, are of particular theo- 
retical interest because they lend themselves to exact so- 
lution in one dimension |^|,^-|^]. The anomalies are typ- 
ically displayed in two ways: one is through the time 
dependence of the reactant concentration c(t) , which for 
the A + A — ► A and the A + A — > reactions decays as 
£-!/ 2 m one dimension instead of the law of mass action 
decay t -1 . The other is through the interparticle distri- 
bution function p(x,t), which is the (conditional) proba- 
bility density for finding the nearest particle at a distance 
x on one side of a given particle. This function scales as 
x/t 1 / 2 , in typical diffusive fashion. In one dimension a 
gap develops around each particle that leads to a more 
ordered spatial distribution than the exponential distri- 
bution implicit in well-stirred systems and "explains" the 
relative slowing down of the reaction. 

In a parallel development, the problem of "anoma- 
lous diffusion" has also attracted a great deal of atten- 
tion ]8|-p0[| . The universally accepted characterization of 
anomalous (as in "not ordinary" ) diffusion is through the 
mean squared displacement of a process x{t) for large t, 
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Ordinary diffusion [a = 1, K\ = D) follows Gaussian 
statistics and Fick's second law leading to linear growth 
of ^:r 2 (i)) with time. Anomalous diffusion is character- 
ized by a nonlinear dependence. If < a < 1 the process 
is subdiffusive or dispersive; if a > 1 it is superdiffu- 
sive. Anomalous diffusion is associated with many phys- 
ical systems and is not due to any single universal cause, 



but it is certainly ubiquitous. Nor is anomalous diffusion 
modeled in a universal way; among the more successful 
approaches to the subdiffusive problem have been con- 
tinuous time random walks with non-Poissonian wait- 
ing time distributions ||, and fractional dynamics ap- 
proaches in which the diffusion equation is replaced by a 
generalized diffusion equation PJfO[|. Some connections 
between these two approaches have recently been clari- 
fied@. 

In this work we consider a combination of these two 
phenomena, namely, the (one-dimensional) kinetics of 
A + A reactions of particles that move subdiffusively We 
pose two questions: (I) How does the reactant concentra- 
tion evolve in time? (2) How does the interparticle distri- 
bution function evolve in space and time? Some aspects 
of this problem have been considered previously using 
the waiting time distribution approach p2 13 1 . The so- 
lutions require approximations relating the reactant con- 
centration to the distinct number of sites visited by a 
particle ]T^| , or the waiting time distributions for sin- 
gle particles to the waiting time distributions for relative 
motion [jl3|| . Here we adapt the fractional dynamics ap- 
proach to the problem and take advantage of the fact 
that the resulting generalized diffusion equations can be 
solved in closed form. We consider both coagulation and 
annihilation reactions. 

Consider first the coagulation reaction A + A — > A 
when the particles move by ordinary diffusion. The prob- 
ability distribution function for the position y of any one 
A particle in the absence of reaction obeys the diffusion 
equation 



d_ 

at 



P(y,t) 
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The coagulation problem can be formulated in terms 
of the probability E(x,t) that an interval of length x 
is empty of particles at time t. This "empty interval" 
function for a diffusion- limited reaction (i.e., one with 
immediate reaction upon encounter) obeys the diffusion 
equation Q 



(3) 



The derivation of this equation is straightforward and 
recognizes that an empty interval is shortened or length- 
ened by movement of particles in and out at either end 
according to the dynamics described by Eq. (g). The 
empty interval dynamics is thus essentially the same 
as that of individual particles in the absence of reac- 
tion, but with a diffusion coefficient 2D that reflects the 
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fact that the relative motion of two diffusive particles 
involves twice the diffusion coefficient of each particle 
alone. The coalescence reaction implies the boundary 
condition E(Q,t) = 1, and E(po,t) — as long as the 
concentration is nonvanishing. 

From E(x,t) one obtains the concentration of particles 



c(t) = - 



dE(x,t) 



dx 



and the interparticle distribution function 

d 2 E(x,t) 
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c(t) dx 2 



(4) 



(5) 



Equation (Q) and the boundary conditions can readily 
be generalized in a number of ways pH , in particular 
to reversible coagulation (A+A^A) and nucleate 
and even to processes with three-site interactions 

The motion of a subdiffusive particle (in the absence 
of reaction) is described by the fractional diffusion equa- 
tion EHil7[ M 
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where o Dt~° is the Riemann-Liouville operator: 
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and K a is the generalized diffusion coefficient that ap- 
pears in Eq. (|lj). Some limitations of this description 
have been discussed recently |flj . 

The construction of the kinetic equation for E(x,t) 
for subdiffusive particles proceeds along arguments anal- 
ogous to those used in the diffusive case. Again, one 
follows the motion of the particles in and out of the ends 
of the empty interval according to the dynamics (||) . This 
readily leads to the fractional diffusion evolution equation 
for the empty intervals for subdiffusion-limitcd reactions 



oD\- a 2K a —E{x,t). 



(8) 



The solution of Eq. (p|) with the boundary conditions 
E(0,t) — 1, and E(oo, tj = can be expressed in terms 
of the Fox H- functions ITq] . In Laplace transform space 
(indicated by a tilde over the function) the solution is 



E(x,u) = - 



dy 
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where s = u a l 2 /y/2K a and E (x) 
Eqs. (Eh and (§) one finds: 
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E(x,0). From 
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where A = c(0) andpo(s) is the spatial Laplace transform 
of the initial interparticle distribution function po(x) — 
p(x,0). 

A commonly considered initial interparticle distribu- 
tion is the random (Poisson) distribution of average con- 
centration A, po(x) = \e~ Xx . For this initial distribution 
c(u) =X/(u+ X^fZKZu 1 - / 2 ) and c(t) is given in closed 
form in terms of the Mittag-Leffler function (h]]21J] of pa- 
rameter a 1 2: 



c{t) = XE a/2 {-\^2K a t 



.a/2 



(11) 



When a = 1 one recovers the usual result for diffusion- 
limited coagulation ^||| since the Mittag-Leffler function 
of parameter 1/2 is E 1 / 2 (—x) = exp(a; 2 )erfc(2;). 

The u —> expansions of Eqs. (g) and ( |i~0| ) are readily 
seen to be independent of the initial distribution po(x) 
and can be Laplace inverted to yield the asymptotic re- 
sults for large t: 



c(t) 



-a/2 



and 



p(x,t)c(t)~ -z H 1 ^ 
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where H is the Fox H-function |9l,ll 
with Eqs. © and (|l|) we find 



'2|. Furthermore, 
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where z = c(t)x is the scaled interparticle distance and 
p a (z)dz = p{x,t)dx. This stationary form is shown in 
Fig. |lj for several values of the diffusion exponent a. 

The interparticle distribution function conveys the in- 
teresting "anomalies" of the problem most clearly. For a 
random distribution of particles on a line this distribution 
is exponential. In particular, the most probable interpar- 
ticle gaps are the smallest. For diffusion-limited reactions 
on a line it is well known that the scaled distribution de- 
viates in two ways from the exponential behavior. First, 
a gap develops around each particle, and the distribu- 
tion vanishes near the origin (see a — 1 curve in the 
figure), indicating an "effective repulsion" of particles. 
Second, the probability of large gaps decays much more 
rapidly than exponentially: the decay goes as a power of 
(9) exp(— z 2 /2). In the subdiffusive case decreasing a leads 
to the diminution of the gap around each particle, that is, 
to a weakening of the effective repulsion and to a behav- 
ior that appears closer to that of a random distribution in 
the short-interparticle-distance behavior. This is evident 
in the progression of the curves with decreasing a shown 
in the figure. Furthermore, the probability of large gaps 



2 



decays as a power of exp(- 



x 



2/(2-o) 



), thus neatly inter- 



polating between the purely random exponential decay 
state as a — > (since then p(x, t) — > c(t)e~ c ^ x ) and 
the more ordered state corresponding to diffusive parti- 
cles at a = 1. Note that the congruence of the curves 
at z = 1/2 is no£ exact although the almost-congruence 
is certainly intriguing (and may be related to the special 
but physically unclear role played by the initial concen- 
tration c eq /2 in the reversible coagulation problem that 
reaches an equilibrium state at concentration c eq , see ||| 
and Eq. @ et seq.). 




FIG. 1. Long-time scaled interparticle distribution func- 
tion for several values of the anomalous diffusion exponent. 
Proceeding upward from lowest to highest curves along the 
y axis intersection: a = 1,0.95,0.8,0.5,0.2. Note that the 
distribution for a — 0.2 on this scale is nearly indistinguish- 
able from the completely random distribution exp( — z) (dot- 
ted curve). 

The arguments presented so far are applicable to the 
irreversible coagulation reaction A + A — > A. The empty 
interval method can not be applied to the annihilation 
reaction A + A — > because annihilation leads to a dis- 
continuous growth of empty intervals. However, recently 
a new method of odd/even intervals has been introduced 
that leads to exact solution in the diffusion-limited case. 
It is based on the construction of an equation for r(x, t), 
the probability that an arbitrary interval of length x con- 
tains an even number of particles at time t [^Q- Again, 
because r(x, t) changes only by the movement of par- 
ticles in or out of the ends of the interval, arguments 
similar to those that lead to Eq. (^) lead to exactly the 
same equation for r(x, t) but with the boundary condi- 
tions r(0,i) = 1 (as for E(0,t)) and r(oo, i)=l/2. The 
concentration of particles is related to r(x,t) precisely as 
in Eq. (@). 

The method of odd/even intervals can again be directly 
extended to the subdiffusive problem, where r(x, t) satis- 



fies the same fractional diffusion equation as E(x, t) with 
appropriately modified boundary conditions. Closed so- 
lution is again possible and, for a random initial distribu- 
tion with average initial concentration A leads to a form 
slightly modified from the result (O): 



c(t) = XE a/2 (-2\^2K~ a t a 
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(15) 



It is appropriate to make contact with the ingenious 
work of Spouge [^J. He introduced a single formalism 
to calculate the particle concentration c(t) for both co- 
agulation and annihilation. His approach is based on 
the probability a(x,t) that two particles, one starting at 
a distance x from the other, have met by time t. In 
ordinary diffusion this quantity also obeys the diffusion 
equation with diffusion coefficient 2D with appropriate 
boundary conditions. Again, his method can be gener- 
alized to anomalous diffusion and the resulting a(x, t) 
is M 
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Spouge's prescription applied to this subdiffussive situ- 
ation leads, for an initially random distribution, to the 
results (|ll]) and ( |l5| ) for coagulation and annihilation re- 
spectively 

The irreversible processes can be generalized in a 
variety of ways for which the empty interval and 
the odd/even interval methods have been suitably ex- 
tended |g|,|U^,|l^] . In particular, consider the generaliza- 
tion to reversible coagulation, A + A # A O. For the 
ordinary diffusion-limited process this generalization is 
accomplished by subtracting from the right hand side of 
Eq. (||) a term proportional to dE(x, t)/dx that describes 
the rate of decrease of E(x, t) due to the fact that par- 
ticles adjacent to either end of an empty interval may 
give birth (at rate v) to a particle that moves into and 
therefore shrinks the interval. The same generalization 
is possible in the subdiffusive case. Assuming that the 
particles born at a rate v move into an empty interval 
subdiffusively, the evolution equation for E(x,t) is given 
by 
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E{x,t) 



d 2 



v-^E(x,t) 



(17) 



This equation can be solved by separation of vari- 
ables jl0|,[l9]]. The diffusive modes have the same spa- 
tial dependence <p n {x) as in the diffusive case, but their 
temporal evolution is no longer exponential. Instead, it 
is determined by the Mittag-Leffler function: E(x,t) = 
'^2ip n (x)E a (—X n t a ) 1 where the A n are the eigenvalues 
associated with the tp n (x). The equilibrium solution 
is E eq (x,t) = exp(-vx/2K a ), 
c eq = v/2K a and p eq (x) = c, 



which in turn leads to 
eq exp(— c eq x), as in the 
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case of ordinary diffusion p]l 4|]. However, the approach 
to equilibrium is algebraic rather than exponential and 
thus qualitatively different, c{t) — c eq ~ t~ a . An inter- 
esting attendant observation involves the dependence of 
these results on the initial concentration c(0). In the 
diffusion-limited case the time constant r in the expo- 
nential decay to equilibrium exp(— t/r) is a function of 
c(0) for c(0) < c e? /2 and changes to a value independent 
of c(0) for c(0) > c e? /2 |. In the subdiffusive case it 
is the prefactor of t~ a that undergoes exactly the same 
change. 

We have considered the coagulation dynamics A+ A — > 
A and A + A ^ A and the annihilation dynamics 
A + A — > for particles moving subdiffusively in one di- 
mension. This scenario combines the "anomalous kinet- 
ics" and "anomalous diffusion" problems, each of which 
leads to interesting dynamics separately and to even more 
interesting dynamics in combination. The fractional dif- 
fusion equation plays a central role in our analysis and 
allows the exact calculation of the density c(t) and of the 
interparticle distribution function p(x, t) within this for- 
mulation. Anomalous diffusion is characterized by the 
exponent a introduced in Eq. (|l|) , ordinary diffusion cor- 
responding to a = 1 . Deviations from ordinary diffusion 
lead to a curious interplay. On the one hand, with de- 
creasing a (and hence increasingly subdiffusive motion) 
the decay of the particle density towards extinction or 
towards equilibrium becomes increasingly slower and in 
this sense increasingly different from law of mass action 
behavior. On the other hand, the spatial distribution 
of initially randomly distributed reactants remains more 
Poissonian for all time as a decreases; indeed, as a de- 
viates from unity the relatively empty regions around 
each particle that are tantamount to (and indeed ex- 
plain) anomalous kinetics in the usual diffusion-limited 
case become more populated. 

A number of generalizations and further questions that 
have been considered in the context of diffusion-limited 
reactions but not for subdiffusion-limited reactions im- 
mediately come to mind. They include an analysis of 
different initial distributions, more detailed consideration 
of reactions with external input (sources), effective ki- 
netic equations for the density, reactions in statistically 
inhomogeneous media, and reactions in the presence of 
external potentials. 
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